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OPTIMAL TIME DECAY RATE FOR THE COMPRESSIBLE 
VISCOELASTIC EQUATIONS IN CRITICAL SPACES 


JUNXIONG JIA AND JIGEN PENG 


Abstract. In this paper, we are concerned with the convergence rates of 
the global strong solution to constant equilibrium state for the compressible 
viscoelastic fluids in the whole space. We combine both analysis about Green’s 
matrix method and energy estimate method to get optimal time decay rate in 
critical Besov space framework. Our result imply the optimal L 2 -time decay 
rate and only need the initial data to be small in critical Besov space which 
have very low regularity compared with traditional Sobolev space. 


1. Introduction and main results 

Many fluids do not satisfy Newtonian law. A viscoelastic fluid of the Oldroyd 
type is one of the classical non-Newtonian fluids which exhibits elastic behavior, 
such as memory effects. The elastic properties of the fluid are described by associ¬ 
ating the fluid motions with an energy functional of deformation tensor U. Let us 
assume the elastic energy is W(U), then the compressible viscoelastic system can 
be written as 

d t p + div(pw) = 0, 

dt(pv) + di v(pv Cg) v) + VP(p) = div(2pV(v) + V(Adiv(u))) 

( L1) dh r ( Wu(U)U T \ 

+Q1V V det(U) ) ’ 

d t U + u-VU = \7uU. 

Here p is the density and v(x,t) is the velocity of the fluid. The pressure P{p) is 
a given state equation with P'(p) > 0 for any p and V{v) = ^(Vu + Ve T ) is the 
strain tensor. The Lame coefficient p and A are assumed to satisfy 

(1.2) p > 0 and A + 2p > 0. 

Such a condition ensures ellipticity for the operator —X7(2pT>-) — V(AV-) and is 
satisfied in the physical case, where A + 2 p/N ss 0. Moreover, Wu{U ) is the Piola- 
Kirchhoff tensor and the Cauchy-Green tensor, respectively. For a special 

case of the Hookean linear elasticity, W(U) = \U\ 2 . 

For the incompressible viscoelastic fluids, there are many important works re¬ 
cently. In [I], the author proved the well posedness problem and find the relation 

X7 k F ij - \7jF ik = F lj ViF ik - F lk ViF ij , 

with F = U — I. This relation indicates that the linear term V x F is actually a 
higher order term. F. Lin, C. Liu and P. Zhang 0|3] proved the local well posedness 
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in Hilbert space H s , and global well posedness with small initial data. In the proof 
of the global part, they capture the damping mechanism on F through very subtle 
energy estimates. At last, in ;|j, the author proved the global well posedness of the 
incompressible version of system CD in the critical L p framework which allows us 
to construct the unique global solution for highly oscillating initial velocity. 

For compressible viscoelastic fluids, in 0 El the authors proved the local and 
global well-posedness in the L 2 based critical Besov type space. Their work used 
the properties of the viscoelastic fluids deeply and their results indicated that the 
deformation tensor U plays a similar role as the density p. It should be mentioned 
that the global existence of a smooth solutions is still an open problem, even in 
for incompressible viscoelastic fluids. P. Lions and N. Masmoudi 7j proved the 
global existence of a weak solution with general initial data in the case that the 
contribution of the strain rate in the constitutive equation is neglected. Recently, 
in [8] , the author proved the global well-posedness in L p based critical Besov space. 

Besides well-posedness theory, optimal time decay rate problem is another im¬ 
portant subject. There are many papers studied optimal time decay rate for com¬ 
pressible Navier-Stokes system 

ranuminHiEiniini. However, due to the 
complexity of the compressible viscoelastic equations, there are little results about 
viscoelastic system. Recently, X. Hu and G. Wu in m give a detailed analysis 
about time-decay rate in the Sobolev space framework. They split the whole sys¬ 
tem into two small systems and then the analysis becomes possible. In [8], the 
author used estimates in homogeneous space and negative Besov space to give a 
slow decay rate when the initial data only small in Besov space with low regularity. 
The main goal of this paper is to get the optimal time decay rate when the initial 
data just small in critical Besov space framework. Hence, we can link the results 
in m and 0 to give a more elaborate characterization about time decay rate for 
compressible viscoelastic system. 

In paper they proved the following proposition which reveal some intrinsic 

properties about compressible viscoelastic equations. 

Proposition 1.1. The density p and deformation tensor U in m satisfy the 
following relations: 



(1.3) \detl7 / 

and U lk ViU ij - U lj ViU ik = 0 

if the initial data (p,U )| t= o = ( po,Uo ) satisfies 



(1.4) 


and U lk ' V, U l 0 J - U l 0 j V ; U ik = 0 


respectively. 

Using Proposition 11.11 the last term in the second equation of m can be 
rewritten as 


(1.5) 
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As in |6j, without loss of generality, we consider Hookean linear elasticity, W{U) = 
\U\ 2 in the following part of this paper. Note that this does not reduce the essential 
difficulties. All results can be easily generalized to the case of more general elastic 
energy functionals. In view of m, we will consider the following system 

{ d t p + div(pv) = 0, 

pdtf + pv ■ W - div(2pV(v)) - V(Adivu) + VP{p) = pW k VjU ik , 
d t U + v • VC/ = VvU, 

(p,v,U)\ t =o = (po,vo,U 0 ) 

where the initial data satisfies m- 

We now state our main result of this paper which gives the optimal L 2 -time 
decay rate for strong solutions in critical Besov spaces. 

Theorem 1.2. Assume that dimension n = 3, p be a constant and I stands for 
the identity vector (1,1,1). There exists 6 > 0 such that if vo € B^{ 2 1 D 
P 0 -P& B 2 / 1 2 n B° iOOI U -I € B n B° hoo and 

II(Po ~ PiU — -OIIb^obJ^ + ^ 

then problem im has a unique global solution {p — p,v,U — I) € C(R + ; B^f) x 

^C(R + ; 1?2( 2 1 ) n L 1 (R + ; ( 2+1 )) x B %( 2 )^ ■ Furthermore, there ex¬ 

ists constant Co > 0, and we have 

(i.7) \\{p — p,v,u — i)(t )\\ B »/’-i < c 0 (i + t)~ n/4 , 

^2,1 

for t > 0. 

Remark 1.3. From we know that the optimal L 2 -time decay rate for compress¬ 
ible viscoelastic equations is 

(1.8) ||(p-p,u,C/-J)(t)|U 2 <C(l+f)-"/ 4 . 

Due to Uj( 2_1 C A 2 , the convergence rate of (E3 is optimal. 

To prove Theorem 11.21 we split the system by Littlewood-Paley operator to low 
frequency part and high frequency part. For the low frequency part, we decompose 
the system into three small system and analyze the green’s matrix carefully as in 
for each small system. Due to the fine properties of homogeneous space and 
singular operators, we can then combine the estimates for small systems together 
to finally obtain an estimates about the whole system. For the high frequency 
part, we reformulate the system as in [5j and using energy estimates in Besov space 
framework to get an appropriate a prior estimates. 

The paper is organized as follows. In Section 2, we introduce the notation, some 
properties of Besov space and some important Lemmas. In Section 3, we split the 
system into three small system and give the estimates for low frequency part. In 
Section 4, we transform the system into an equivalent form and prove an estimates 
for high frequency part. In Section 5, we give the proof of Theorem 11.21 
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2. Preliminaries 

In this section we first introduce the notation which will be used throughout this 
paper. Secondly, we give some basic knowledge about Besov space. At last, we 
present some useful Lemmas and Theorems. 

2.1. Notation. Let n stands for the dimension, L p (l < p < oo) denote the usual 
L p -Lebesgue space on 1". [ 2 ] stands for the integer part of a number 2 G R. The 
inner-product of L 2 is denoted by (•,•)• If S is any nonempty set, sequence space 
£ P (S) denotes the usual l p sequence space on S. For any integer £ > 0, X7 e f denotes 
all of £-ih derivatives of /. 

For a function /, we denote its Fourier transform by IF[f] = f: 

F[M) = m = (2 K)- n/2 f f(x)e~ ix <dx. 

J R n 

The inverse of IF is denoted by J r_1 [/] = /: 

^- 1 [/](*) = fix) = (27T)-"/ 2 [ me* x dt. 

JR" 

2.2. Besov spaces. In this section, we will give some basic knowledge about Besov 
space, which can be found in |19] . First we introduce the dyadic partition of 
unity. We can use for instance any (<j>, y) G C°° , such that <j> is supported in 
{£ G r : 3/4 < |£| < 8/3}, X is supported in {(gl": |£| < 4/3} such that 

x(0 + E^ 2_9 0 = 1 ?6R n , 

9>0 

E</(2- 9 o = i if 

Denoting h = J-~ l [(j)\ and h = J r_1 x, we define the dyadic blocks as follows 
A_iit = x(D)u = h* u, 

A q u = cp(2~ q D)u = 2 qn f h(2 q y)u(x — y)dy if q > 0, 

JR” 

A q u = (j){2- q D)u = 2 qn [ h{2 q y)u(x - y)dy if q G Z. 

J R” 

The low-frequency cut-off operator is defined by 

SqU = E Sq u = E ^ kU ■ 

— l</c<q—1 k<q— 1 

The formal two decompositions 

u = A q u, u = E A q u 

q> — 1 gEZ 

are called inhomogeneous and homogeneous Littlewood-Paley decomposition re¬ 
spectively. 

Let us give the definition of inhomogeneous Besov space as follows. 
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Definition 2.1. For s £ 1 and 1 < p, r < oo, and u G S'. The inhomogeneous 
Besov space B 8 r consists of distributions u in S' such that 


Mis* 



< +oo. 


Let us now introduce the homogeneous Besov space. 


Definition 2.2. We denote by S' h the space of tempered distributions u such that 

lim S q u = 0 in S'. 

q — y —oo 


Definition 2.3. Let s be a real number and (j>, r) be in [1, oo] 2 . The homogeneous 
Besov space B 8 r consists of distributions u in S' h such that 

Mbj iP “ (E 2rj 'i A ^n w 

\«ez 


i /r 


< +oo. 


From now on, the notation Bp, B 8 will stand for B 8 ^ and B 8 X respectively. The 
notation B 8 , B s will stand for 1 and x respectively. 

The study of non stationary PDE’s usually requires spaces of type L r T (X) := 
L r (0,T;X) for appropriate Banach spaces A'. In our case, we expect A to be a 
Besov spaces, so that it is natural to localize the equations through Littlewood- 
Paley decomposition. We then get estimates for each dyadic block and perform 
integration in time. However, in doing so, we obtain bounds in spaces which are not 
of type L r (Q,T]Bp) or L r (0,T; B 8 ). This approach was initiated in [18] naturally 
leads to the following definitions for the inhomogeneous Besov space. 

Definition 2.4. Let (r,p) 6 [l,+oo] 2 , T G (0,+oo] and set. We set 

IMIlj,(B‘) := 



and 

Lr(Bp) := e IMI< +°°} ■ 


Owing to Minkowski inequality, we have L!j,(B 8 ) °A Lj,(B 8 ). That embedding 
is strict in general if r > 1. We will denote by Ct(B 8 ) the set of function u 
belonging to L^(B 8 ) 0 C([0, T]; B 8 ). For the homogeneous Besov space, we can 
define similarly. 

Let X stands for B or B , we have the following interpolation inequality: 


Mil 


p(A-) ^ IMIl^X* 1 ) 


L r J(x; 2 y 


with 


- =-1-and s = 9si + (1 — 0 )s 2 , 

r r i r 2 

and the following embeddings 


L r T (X;/ p ) c —>• L r T (C 0 ) and C T (X£ /p ) -a C([0, T] x IT). 
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Another important space is the hybrid Besov space, we give the definitions and 
collect some properties. 

Definition 2.5. let s,t £ R. We set 

■•= E 2 9S ||A ?U || L9 + E 2 qt \\\u\\ LP . 

q<Ro q>Ro 

and 

B s q f p (^ N ) := {u e S' h (R N ) ■■ Mb*.* < + 00 } , 

where Ro is a fixed large enough number determined in the proof of global existence. 

Lemma 2.6. 1) We have B 22 = B s . 

2) If s <t then B^ p = 5®nBj. Otherwise, = B p + I?*. 

3) The space B p ’ p coincide with the usual inhomogeneous Besov space. 

4 ) If s 1 < S2 and ti > t2 then B p Wi B p 2 ^ 2 . 

5) Interpolation: For si, S 2 , <7i, 01 £ M and 8 £ [0,1], we have 

||/|l B ^ 1 + t 1 -e)»2.^ 1 + ( 1 -9)"2 < ll/ll^ 1 ^ 1 11/11^2^2 • 

From now on, the notation B^ will stand for _£?*’* and the notation B s,t will 
stand for B^- For more information about Besov space and hybrid Besov space, 
we give reference P2J 2D] I2T] [25]. 

In the last of this introduction, for the reader’s convenience, we list an important 
Lemma EDEI which will be used in the following. 

Lemma 2.7. Let F be a homogeneous smooth function of degree m. Suppose 
1 — n/2 < p < 1 + n/2 and —l/n < p' < n /2 + 1. Then the following inequialities 
hold: 

\(F(D)A q (vVc)\F(D)A q c)\ 

< Ca q 2~ q ^ p ~ m ^ Ill’ll Bn/ 2+1 ||c|| bp' ||F(I?) A g c||i, 2 , 

\(F{D)A q vVc\F(D)A q c)\ 

< (7a g 2-^- m )min(2«,l)||r ) || i j n/2+1 ||c|| B) ,-i, P ||F(D)A gC || i 2, 

\(F(D)A q (v ■ Vc)|A,d)| + |(A,(v • d)\F(D)A q c)\ 

< Ca 9 2"‘ ?(p " m) min(2 9 ,l)||?;|| i j„ /2+1 (||c|| B p-i,p||A 9 d|| i 2 + ||d|| B p-i, P ||A 9 c|| L 2 ), 

\(F(D)A q (v ■ c)\A q d)\ + |(A,(t, • X7d)\F(D)A q c)\ 

< Ca q \\v\\g n / 2 +i(2~ qp ' ||F(D)A g c||i 2 ||d||^p/ 

+ nrin(2 9 , l)||d|| B ,-i„||A,dM. 

2.3. Useful Theorems. In the part, we will list two Theorems about well-posedness 
of equations (11.61) which are essential for our proof of Theorem 11.21 Denote 

E P (T) := {t; e U([0, T];B p ^ p ), d t v,V 2 v £ L\0, T; 

For v £ E p (T) will be endowed with the norm 

IMIe p (t) := lkll L == (B ^- 1 ) + W dtV ’ v 2 v h} r (B? /p - 1 y 
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Through similar methods used in [8] or just change the low frequency estimates in 
[5], we will obtain Theorem 12.81 Due to the proof has no new ingredients, we omit 
it. 

Theorem 2.8. Let 1 < p < In and n > 2. Let vo be vector field in BfJ v 1 . 
Assume that po satisfies ao := po — 1 £ B p ^ p and Uo satisfies Fq := Uo — I £ B'p V 
and 


(2.1) infp o (ir)>0. 

X 

Then system em has a unique local solution (p, v, t/) with v £ E p (T), U — I £ 
C([0, T]; Bp^ p ), p bounded away from 0 and p—l£ C([0, T\; B p ^ p ). 

Denote: 

£‘ 5 “ {(o,«, F) £{L\ 0, oo; B^ p +1 ' s ) n Z°°(0, oo; B 8 ^ 1 ' 8 )) 

x (L\0, oo; B s 2 * p +1 ’ s+1 ) n Z°°(0, oo; B?" 1 ’'" 1 ))" 
x (T 1 ( 0 , oo; n T°°( 0 , oo; 

where s p = s - ^ + §. 

The global well posedness of equations m in Besov space framework are as 
follows. 


Theorem 2.9. [8) Let p > 0 be a constant such that P'{p) > 0. Suppose that 
n = 3. There exist two positive constants ao and C such that for all (po,Vo,Uo) 

•i 7 - C) n/ 2 — l,n/p T r T — T->n/ 2 —l,n/p _ T->n/2—l,n/p—1 , 

with Po - p e B 2 ; , Uo - I G B 2 f p , Vo G B 2 ' , and 


( 2 . 2 ) 


|| PO — P || R n/2-l,n/p + Il'Poll R n/2-l,n/p-l + 11 C/q ~ /| | R n/2-l,n/p < Op, 


t/ien if 2 < p < n, system EM has a unique global solution (p — p,v,U — I) £ £ n / p 
with 

II (P- P,V, U - I)\\ £ n/p < C(||p 0 - p|| „/2-1,„/ p + \\V 0 \\ „/2-l,n/p-l 

\ 2 ,p 2 ,p 

+ || — U\\gn/2-l,n/p\ . 

^2 ,p / 

Remark 2.10. Taking p = 2 in Theorem 12.91 we will get global well-posedness in 
the critical homogeneous Besov space framework. Assume 


|| (PO P,U ^)Hsn/2 nR 0 ^ + ||'Po|lsn/2- 1 nBO ^ < 5, 

as in our main Theorem ll.2l If S > 0 is taken to be small enough, the above assump¬ 
tion will imply (12.21) . hence, we obtain the results in Theorem 12.91 Particularly, we 
know 

roo 

(2.3) / \\v\\g n / 2 +idt < CS. 

Jo 


This estimate plays an essential role when we estimate the nonlinear terms. 
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3. Analysis about low frequency part 


In this section, we first decompose the system into three small scale system, then 
analyze the semigroup carefully for the low frequency. Without loss of generality, 
we assume P'{ 1) = 1, p = 1 and set v = A + 2/z, A = fiA + (A + ^i)Vdiv. Define 

K(a) = P ^ 1 + a ) - 1, d= A -1 divu, 

1 + a 

Si = A _1 curk> with (curlii)ij = d Xj v l — d Xi v J , 

£ij = A~ 1 d Xi A~ 1 d Xj (F ij + F ji ), 

W = A ~ 1 d Xk (F lk \/iF lj - F l: ShF ik ) - A~ 1 d Xk {F lk S7iF ji - F li S/iF jk ). 

where 

A7 = ^- 1 (l?r/) for s £ R.. 

Performing same procedure as in [|B], we will obtain 

dta + A d = L — v ■ S/a, 
dtd — fiAd — 2A a = G — v ■ S/d, 

(3.1) d t S + 2Ad = J-v- S/S, 
d t (F T — F) + ASl = I — v ■ S7(F T - F), 

d t n - fxAn - A(F T — f) = h — v ■ vn, 

where the equation about d have the following equivalent form 

(3.2) d t d — vAd — AS = K — v ■ S/d, 
where 


L = — adiv?;, 

G = v ■ Vd + A _1 div (—v ■ S7v + FS/F — K(a)S7a - -—Av — div(aA 1 )), 

1 + a 

H = v ■ S7SI + A _1 curl (—v ■ S7v + FS/F — K(a)S7a -—— Av) + W, 

1 + a 

I = ( S/vFf - S/vF, 

J = — [A~ 1 d Xi A~ 1 d Xj , v k ]d Xk (F^ + F*) 

+ A- 1 d Xi A- 1 d Xj {(S/vF) i ^ + ( S/vFY i ), 

K = v ■ S/d + A~ 4 div (—v ■ S7v + FS/F — K(a)S7a - -—Av + div(aF)). 

1 + a 


Here, we denote 


Mr(t) 

M 2 (t) 


sup (1 + r)"/ 4 (||a(r)|| fl „ /a -i,« /a + ||rf(^)||^„ /a _ 1 ), 

0<r<£ 


sup (1 + T) n ' 4 {\\S(r)\\ B n/2-i,n/2 + ||d(r)||^ B/a _ 1 ) , 
0<r<i 


Ms(t) := sup (1 + r)"/ 4 (||(F T - F)(r) || s „ /a _x.„ /a + 11^(711^-0 , 

0<r<£ 


M 4 (t) := sup (1 + r)-/ 4 (||a|| i2 + ||F|| L2 + ||u|| L2 ), 

0<r<£ 
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and 


M{t) sup (1 + t)"/ 4 (||a|| B n/2 + ||F|| B n/2 + IM|b«/ 2 -i) • 

0 <T<t 


From the basic properties of Besov space, we easily know 

M(t) « Mi(t) + M 2 (t) + M 3 (t ) + M 4 (t), 

under the smallness condition of initial data, where we used (5.13) to (5.15) in [8]. 


We here also denote 



(3.3) 

Mi = 

/ L — v ■ 
\G-v 

■ Va 

■Vs 

and 




(3.4) 


m 3 = 

(’ 


J — v ■ VS 
K — v ■ Vd 


M 2 = 


H — v ■VO 

Now, we need to introduce the following linearized system with convection terms. 

<9 t a + Ad = L — u • Vo, 

<9 t d — pAd — 2Aa = G — v ■ Vd, 

(3.5) <9 t £ + 2Ad = J — v ■ VS, 

<9 t (F T — F) + Aid = I - v ■ V{F T - F), 

- /rAfl - A (F T - F) = H-v VQ. 

We can decompose the above system into three subsystems. 


(3.6) 

J d t a + A d = L — v ■ Va, 
yd t d — jiAd — 2A a = G — v ■ Vd. 

(3.7) 

j d t S + 2Ad = J —v ■ VS, 

| dtd — vAd — AS = K — v ■ Vd. 

(3.8) 

j d t (F T -F) + Afl = I-v V(F t - 

1 d t n - pAn - a (f t — f) = h — v ■ 


It is easily observed that the above three system are similar, so we now only 
study the following linear system. 


(3.9) 


d t c + aAu = 0, 
d t u — kAu — f3Ac = 0, 


where c, u are scalar functions and a, /3, n are positive constants. We first give 
some important properties of the Green’s matrix for the above system. 

Lemma 3.1. Let Q be the Green matrix of system m - Then we have the following 
explicit expression of Q: 


s&t) = 




A+-A_ 

A_l £ A_£ 


-A_ 


lei 


—a 

x + 


( A_i_£ A _ £ \ . 

^-) 1^1 


A + -A_ 

A + t_ A _ e A_ I 


A + -A_ 


a± = --k,|ci 2 ±-\AfFW- 


where 
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Lemma 3.2. Given R > 0, there is a positive number 6 depending on R such that, 
for any multi-indices 7 and |£| < R, 


7l 


\Dp(U)\<Ce- e ^ 2t (i + \t;\)^(i + t)\ 
where C = C(R , j'y|). 

The proof of the above two Lemmas follows from Lemma 3.1 and Theorem 3.2 
in [53], so we omit the proof for simplicity. Next, we prove an important Lemma 
which plays the key role in the low frequency analysis. 

Lemma 3.3. Let Q be the Green matrix of system (EM, dimension n = 3. Denote 
U 0 = (c 0 ,u 0 ), then Q(t ) satisfies the estimate 

E \m)\u 0 \\ L * < c(i+ t)- n / A \\u 0 \\go ^ 

q<R 

for t > 0 and R > 0 is a large positive constant. 

Proof. By Placherel’s theorem and Lemma 13.21 we have 


1/2 


wgm.uoh 


< 


< 




df 


( [ e~ e ^ 2t 

MOUo 

y f2«<ici<f29 



1/2 


(3.10) 


1/2 


<||A 9 I7o|Ux / 

|28<|g|<|2« 

<t-"/ 4 ||A,^olUi ' f 


e-^df 


r 2 e 


-Br‘ 


dr J 


1/2 


V*j29<|4|<Vt§2« J 

1/2 


<t~ n/A \\U 0 \\^ lao e ~ iAHe 


1 j 2 

Now we do some calculations to bound J2 q <R e ~^ 4,lte (/ ~ . Let k 

—q, then we have 


^ e -!«(i- e -f w ) 1/2 = E e~^T te (l-e-f (i)** 6 ) 


1/2 


q<R 


k——R 

=1 +11 + HI, 


where 


1= E e-i(^) fct0 (l-e-f U)^) 1/2 , 

k=-R 

[log 4 


11= E e~i^) kte (l-e~^) kte ) 

k =1 


1/2 
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and 


111= jr e-^) kte (l-e-™(iT t0 ) 1/2 . 

k=[iog 4 ^te]+i 

For I, we have I < R • 1 < C < oo. Due to for arbitrary t > 0, there exist positive 
integer N > 0 such that tO < ^4^. Without loss of generality, we can choose 
tO = ^4^. For 77, we have 

N 


Z/=£e- 

k=1 

AT AT 

<C^e-^(3) fe4JV <C'^e-s 


_ l_A-(k-N) 

30 ^ 


k=l 

N—l 

<CJ2 e~^ 4m < C 

m =0 

By Taylor’s formula, we have 

„ 20 pVjg 20 /l\ fe „ ^ 

i_ e - T w f»+E 


fe=l 


< OO. 


(IQ 


n! 


When /c > log 4 ppfd, we have ^ (j) td < 1. So by the properties of alternating 
series, we know 


(3.H) 


20 1 1 ISa 40 /1 \ 

1 _ e --U,) tB < — - tO. 

~ 3 V 4 / 


Substituting (13.111) into III , we obtain 


fe=l 


hi < c y e ~ mte (j) ^ 


4 / 


[iog 4 fte] k 

<CVtOe~^ te Y 

k=i 

<C < oo. 

Combining the estimates for I, II, III, we finally get 
(3.12) ' ' V2 


^ e -i 49te (l^e-^ e ) < C<oo, 


q<R 

where C does not depend on t. 

Summing up estimate (13.101) with (13.121) . we obtain 

E WGmML* <Ct- n ^\\Uo\\so ioo E«" |4,fe (l^e-f 4<,te ) 
(3.13) q<R q<R 

<< 7 t- n/ 4 ||£/o|U . 


1/2 
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Similarly, we also find that 


j2\\Gm q u 0 \\ L z <cj2\\\u 0 \\ L 

(3.14) i<R <i<R 

(Vsy < C < oo. 


n—l—6r z t 

r e 


|29<r<§29 


\ 1/2 
dr J 


q<R 


Combining (13.131) and (13.14[) . we finally arrive our desired results. 


□ 


Remark 3.4. Denote Q\, C/ 2 and G 3 represent the Green matrix of system (EH), EH 
and (13.81) respectively. Denote Vq = ( ao , do), Vq = (£ 0 , do) and V§ = {Fq—F 0 , D 0 ), 
then using similar methods as in the proof of Lemma 13.31 we will have 

£ < ca + tr^^Wzo^ 

q<R 

for i = 1 , 2 ,3. 

Next, we need to consider estimates about M 4 , M 2 , M 3 defined in EH and 

EH- 


Lemma 3.5. There exist an 8 > 0 such that if 

ll a o||,B'*/2 + ||^o||s"/2-l + ||Fo|| B n/2 < 6 , 

then there exists a constant C > 0 independent of time T such that 

\\M\, M 2 , M 3 ||^o ^ < C(l + t)- n / 4 M{t)f{t)+C(l+t)~ n ' 2 M 2 {t ) 

fort e [0, T\, where f{t) = \\v(t)\\g n/2 +i e L^O, 00 ). 


Proof. Now, we start with M 4 . For term v • Va, we have 

||v - Va||^o_ <C||w Vo||£i < C\\v\\ L 2 \\Va \\ L 2 


(3.15) 


< C (1 + t)- n/ 4 M 4 (t)||a|| B n/ 2 -l,n /2 
<C(1 + t)~ n/2 M 2 (t). 


For term adivu, we have 


ll a divu|| B o ^ <C'||adivi;|| L i < C \\a\\ L 2 1| Vv || L 2 


(3.16) 


<C(l+t)~ n/4 M 4 (t) £ || A , Vt ;|| L 2 + ^ || A , Vv ||£2 

\q<R q>R 

<(7(1+ i)-” / 4 M 4 (i) (||u || L 2 + |H| Bn/2+I ) 

< C (1 + t)~ n ' 2 M 4 (t) (M 4 (i) + f(t)) 

< C (1 + t)~ n/2 M 2 (t) + C (1 + t)~ n/4 M(t)f(t). 


For term v • V d, we have 

\\v- Vdllso^ <C7||w Vd||ii < C||u|| L 2 ||Vd || L 2 

<C\\v \\ L 2 (||v|| L 2 + ||u|| B „/2+l) 

< C (1 + t)~ n/2 M 2 (t) + C (1 + t)- n,i M{t)f(t). 


(3.17) 
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For term A 1 div(i; • Vu), we have 


(3.18) 


||A 1 div(i; • Vu)||< C ||u • Vu||^o < C ||u • Vu|| L i 

<(7(1+ t)~ n/2 M 2 (t) + C (1 + t)~ n/4 M(t)f(t), 


where we used similar argument in (13.171) to get the last inequality. For term 
A _ 1 div(FVF), using similar method as in (13.181) and (13.151) . we have 

(3.19) ||A- 1 div(FVF)|| bo 1oo < C{1 + t)~ n/2 M 2 {t). 

Using composition rules (for example: Theorem 2.61 in (19) 1 and similar argument 
as above, we obtain 


(3.20) 


A _1 div 


1 + CL 


Av 




<C(1 +t)~ n / 2 M 2 {t) 

+ C{l + t)~ n/4 M(t)f(t). 


Summing up estimates (13.151) to (13.201) . we get 

(3.21) ||Mi||s ?oo < C(l+t)- n/2 M 2 {t) + C(l+t)~ n / 4 M{t)f(t). 


Next, let us due with Fi. The term v- V£, v-A/d and I\ all can be estimated similar 
to the term appeared in Fj, so we just need to give the following estimates about 
J. Since 


(3.22) 


A-^A^V*. [(A7vF) ij + (VvF) ji \ ^ q < C'||Vu|| i 2 ||F || i2 
< C (1 + t)~ n/2 M 2 {t) + C (1 + t)~ n/4 M{t)f(t), 


|| \A~ l d Xi AT 4 d X;j ,v k ] d Xk (F^+F^\\ <ao 
<||A^ i A- 1 ^(^5 Xfc (F^+F^))|| J j U 

+ \\v k ■ A~ 1 d Xi A~ 1 d Xj (d Xk (F ij + ^ i ))||^„ ^ 
<C(l+t)- n ' 2 M 2 (t) + C\\v\\ L 2\\A- 1 d Xi A- 1 d Xj (d Xk (F i i+Fi i ))\\^ 
< C (1 + t)~ n / 2 M 2 {t) + <7 ||w|U» (||-F||l= + ||F|| B „ /a -i,» /a ) 

<C{l + t)~ n/2 M 2 (t), 


we have 


(3.24) W M 2\\b° 1oo < C (l + t)- n ' 2 M 2 (t) + C {l + t)~ n / 4 M(t)f(t). 

Due to all terms appeared in M 3 can be estimated similar to the terms appeared 
in Mi and M 2 , here, we just give the estimates as follows 

(3.25) IIMall^ < C(l + f)-"/ 2 M 2 (f)+C(l+t)-"/ 4 M(t)/(t). 

At this stage, we easily finished the proof by just summing up (13.211) . (13.241) and 
(1+251) . ' ' ' □ 


Denote V = ( a,v,F ), V\ = (a,d), V2 = (£, d), V3 = ( F T — F, ft) and define 
A jif '■= A_i/ + Yha<q<R^qf f° r a tempered distribution /. Now, we can prove 
the following proposition which is the main results for low frequency part. 





























14 


J.X.JIA AND J. PENG 


Proposition 3.6. Let n = 3, there exists an <5 > 0 such that if 

||a 0 || B n /2 + HuoIIb "-/ 2 - 1 + ll-^olls "/ 2 < ^ 

then there exists a constant C > 0 independent of time T such that 

sup (l+r) n ^\\A R V(r )\\ L 2 < C\\V 0 \\go + CSM(t) + CM 2 (t) 

0 <T<t 1 ’°° 

for t G [0, T\. 


Proof. By the properties of Littlewood-Paley operator and (5.13), (5.14), (5.15) in 
[8], we have 

\\A r V(t)\\ l2 < \\A q A R V(r)\\ L2 < \\\V(t)\\ l * 

q<R 

(3 26) ~ E II AgVi ( r ) II A 2 + E II AgV^ (t) II L 2 

q<R q<R 

+ J2 l|A g P 3 (r)|U 2 + E ||A 9 A- 1 (FVF)(r)|| i2 

q<R q<R 

For the last term appeared in the above inequality (13.261) . we have 


(3.27) 


E ||A 9 A- 1 (FVF)(r)|| i 2 <C\\F\\ 6 o 2 \\F\\ 6 n/, 

q<R 

<C(1 + t)~ n/2 M 2 (t), 


where we used Lemma A.4 in | 8 ](take t = s = 0, s = t = ^,p = 2 and 7 = 0). 
From EH and E3, we easily get 

AgVi(t) = Qi(t)A q Vio + j Gi(t — s)A q Mi(s)ds, 

Jo 

(3.28) A q V 2 (t) = g 2 (t)A q V 20 + I G 2 (t - s)A q M 2 (s)ds, 

Jo 

A q V 3 (t) = g 3 (t) AqV 30 + [ g 3 (t - s)A q M 3 (s)ds. 

Jo 

By using Remark 13.41 and Lemma 13.51 we will get 

W^rV(t)\\ L 2 < J 2 HA 9 Vl(r)|| i 2 + £ ||A g P 2 (r)|| L 2 

q<R q<R 

(3-29) + £ ||A g V 3 (r)||i 2 + (1 + r)-"/ 2 M 2 (r) 

q<R 

</ + //+( 1 + T)~ n / 2 M 2 {r), 

where 

I = E {ll0i(T)A 9 R 1 (T)|| i 2 + \\g 2 {r)A q V 2 (T )\\ L 2 + ||e 3 (r)A^3(r)|U a } 

q<R 

A (1 + t) - "/ 4 {||Vio||sj ^ + H^2 0 |lsJ jOO + ll^°Hsf i00 } 

<C (1 + t)-"/ 4 ||Ro|Ibo , 

l,oo 


(3.30) 
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II = E r{mr - *)A,Vi (*)|| i3 + II g 2 (r - s)A,V 2 ( a )|| L 2 

q<R^° 

+ \\@3(j~ - s)A g F 3 (s)|| i 2|ds 

< C j\l + r - S )- n/4 {|| VlII^o _ + IIVSHbo _ + IIVsIIbo^} ds 
<C(l+T)~ n/ 4 M(t) [ f(s)ds + C (1 + T )~ n/2 M 2 (t) 

Jo 

< C (1 + T)- n/4 M{t)S + C (1 + r)" n/2 M 2 (r). 


Combining (13.291) . (13.301) and (I3.31|) . we finally obtain 


(3.32) 


sup (l + r) n / 4 ||AflC(r )|| L 2 < C||C 0 |Ibo + C5M(t) + CM 2 (t). 

0 <T<t 1 ’°° 


□ 


4. Analysis about high frequency part 

In this part, we need to transform the equation into another form and estimate 
in the high frequency domain which is completely different to the method used for 
the low frequency domain. 

Without loss of generality, assume p = 1 and 7 = i/P'(p) — 1. Denote a = p — 1, 
F = U - I, A = (—A) 1 / 2 , d = A~ 1 divi;, e ij = From U ek \/gU ij - 

t/^V gU lk = 0, we easily know 

A _1 (VjVfcP* fc ) = -A F ij - A -1 Vfc (F ej VgF ik - F ek V e F ij ) . 


Hence, we can transform the equations (ESP into the following new form. 
d t a + v ■ Va + Ad = G\, 
d t e ij + v ■ Ve ij - pAe ij - (A + p)V t Vjd 
1) + A -1 ViVya + A F ij = Gl j 

d t F ij + v ■ V/ ' J - Ae ij = Gl j , 

where 


and 


Gi = a div-y, Gl j =V k v i F kj , 


G% =v ■ Ve ij - A _1 Vj 


v ■ Vv i + C(cl)Av + F'' Y , F' k 


+ A - 4 V fc (F^ VF ifc - F ek VgF ij ) 


with C(a) = K(a) = P '[]^ a a) - 1. Moreover, we have 


(4.2) = -Vj-o + Gl G{ = -Vi(aF^). 

Now, we give the main estimates for high frequency domain in the following 
proposition. 
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Proposition 4.1. There exists an 6 > 0 such that if 

ll a o||s™/2 + Ikolls"/ 2 - 1 + II-FoIIb"/ 2 < A 


then there holds 

j t E q {t) + c 0 E q (t) < C {a,(l + i)~ n ' A M(t)f(t) 

+ ot q \\Gi, G3|| B n/2-i,r>/2 + a g ||Go, G 2 ||b-/2-1 | 
for t G [0, T] and q > R, where a q < 1, 

/*oo /*oo 

/ f(t)dt= \\ v (t)\\g n/2+1 dt < CS 
Jo Jo 


and co dose not depend on q. Here, E q (t) is equivalent to 2 2 3|| A (? a||_£ / 2+2 2 9||A g F||^2 4- 
2(f- 1 )9||A g e|| B 2 . That is, there exists a D\ such that 

——E q < E q < D\E q 
D 1 

where 

E q = 2^\\A g a, A q F \\ L 2 + 2^~ l )<i\\A q e\\ L 2 


Proof. Applying the operator A q to system (EH), we find that (a, e, F) satisfies 
A q d t a + AA q d = A q G\ — A q (v ■ Vo), 

A q d t e ij - /. iAA q e ij - (A + /x)ViV 7 -A„d + A _1 V i V.A„a 
(4.3) ... 

+ AA q F lJ = A q G% -A q {v Ve lJ ), 

A g d t i™ - AA,e« = A ? g^' - A,(u • vf«), 

where i,j = 1,2,3. Taking the L 2 -product of the second equation of (14.31) with 
A q e 13 , then summing up the resulting equation with respect to indexes i,j , we can 
get 

(4 4) 2 ~dt + l 2 + (^ + M)l|AA ? d|| 2 2 - (A g a|AA 9 d) 

+ (AA q F\A q a) = (A q G 2 \A q e) - (A q (v ■ Ve)|A,e), 

where we used the fact d = —A _2 VjVje*^. We apply the operator A to the first 
equation of (14.31) and take the L 2 -product of the resulting equation with —A q d, 
and take the L 2 -product of the second equation of (14.31) with A -1 V;V_, A g a. Then, 
summing up the resulting equations yields that 

- j t (AA q a\A q d) - ||AA,d|| 2 2 + ||AA,o ||£2 - (A + 2 M )(A 2 A 9 d|AA 9 o) 

( 4 - 5 ) + (AgF^IVjVjAgo) = -(AA g Gi|A 9 d) + {AqG^A-^y.Aqa) 

+ (AA q (v ■ Va)\A q d) - (A q (v ■ Ve y )|A -1 ViVjA,a). 

We apply the operator A to the third equation of (14.31) and take the L 2 -product of 
the resulting equation with A q e^ and take the L 2 -product of the second equation 
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of (14.31) with A A q F % F Then, summing up the resulting equations yields that 

^(AA q F\A q e) - ||AA 9 e|| 2 i2 + ||AA 9 F || 2 2 + /r(A 2 A g e|AA 9 F) 

(4 6) + (A + n)iy i V j A q d\AA q F ij ) + {V i 'V j A. q a\A q F ij ) 

= (A q G 2 \AA q F) + (AA q G 3 \A q e) - (AA q (v • Ve)|A q F) 
~(AA q (vVF)\A q e). 

Now, applying the operator A to the first and the third equations of (14.31) . then 
taking the L 2 product of the resulting equations with AA q a and AA q F 1 ^ we will 

i|||AA g a|| 2 2 + (A 2 A 9 d|AA 9 a) 

= (AA g C?i|AA 9 a) — ( AA q (v • Va)|AA 9 a), 

i|||AA,F|| 2 2 -(A 2 A g e|AA ? F) 

= (AA,G 3 |AA g F) - (A,(t; • VF)|AA 9 F). 

We apply the operator A _1 V;Vj to the third equation of (14.31) and take the sum¬ 
mation with respect to i,j, then we take the L 2 times the resulting equation with 
A _1 ViVjA g i™ to get 

i|||A- 1 V i V,A 9 F ^||| 2 + (AA q d\ViVjA q F ij ) 

(4 - 9) = (A- 1 V l V i A 9 G^|A- 1 V fc V,A q F k ’ e ) 

- (A-^iVjA,^ • VF ij ')|A -1 V fc ViAgF M ). 


get 

(4.7) 

and 

(4.8) 


Summing up (14.41) . (14.51) . (14.61) and (14 . 7[) - (IT.9D yields that 

\j t f 2 q +f q + 2(A g F ij | VjVjAgo) = (A q G 2 \A q e) - (AA q Gi\A q d) 

( 4 10) - (A 9 |A -1 ViVjA 9 a) + (AqG 2 |AA 9 F) + (AA 9 G 3 |A 9 e) 

+ (A + 2/i)(AA 9 Gi|AAqa) + /j,(AA q G 3 \AA q F) 

+ (A + ^(A-^iVjAgGalA-^fcV/A,^) + F q , 

where 

fq = ll^<? e lll 2 + (A + 2/x)||AA g a||| 2 + M||AA 9 F||| 2 

+ (A + /x)||A- 1 V i V J -AqF^|| 2 2 - 2{AA q a\A q d) + 2(AA q F\A q e), 

fq = (f‘ — l)l|AA 9 e||| 2 + (A + n — l)||AA g d|| 2 2 + ||AA g a||| 2 + ||AA 9 F||| 2 
- (A 9 a|AA 9 d) + (AA 9 F|A 9 e), 
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F q =- {A q (v ■ Ve)| A g e) + ((AA,(v • Va)|Agd) + (A,(« • Ve^A-WiV^a)) 
+ ((AA,(v • Ve)| A g F) - {AA q (v ■ VF)|A,e)) - m(AA,(« ■ VF)|AAgF) 

- (A + ^XA-^VgAg^ • VF)|A- 1 V fe V^AgF M ) 

- (A + 2fi) 2 (AA q (v ■ Va)|VAga) 


Here, we can take R to be a fix large enough constant. For q > R, we can easily 
deduce 

(4.11) f q « 2 29 || AgO||^2 + ||Age||| 2 + 2^||AgF||| 2 , 
and 

(4.12) 2 29 ||A g a || 2 2 + ||A,e||£ a + 2 2 ^||AgF|| 2 2 < ~f q . 

Using the identity (14.21) . we find that 

(A,F” V,V,«) =(V l V i A q F i XAg a) 

(4.13) 

= ||AAga||i2 + (AAga|A- 1 V J G J 0 ). 

Let E q {t) = 2^~ 1 ') q f q , then we have 

(4.14) Eq (t) » 25«||Aga|| i2 + 2 (^- 1 )^||A q e|U 2 + 2*«||A,F|| ia 
By (14.101) . (14. 1 1I) - (I4.14I) and Lemma [2771 we finally obtain 

j t E q {t) + c 0 E q (t) < Ca q (l + + CaollGx, G 3 || B „,2-i,„,2 

+ Ga,||G 0 ,G2||B»/ 2 -i. 


□ 


5. Derive optimal time decay rate 


With the analysis about low and high frequency part, we now give the proof of 
Theorem P From Proposition 14. 1[ we know that 

r t 


(5.1) 


E q (t) <e- Cot E q ( 0 ) + C e- c °^(a q {l + t)-"/ 4 M(t)/(t) 

+ <Tg||Gl, G 3 || B n/2-l,n/2 + CIg||Go, G 2 11 B n/ 2-1 ^ d"?". 

Through homogeneous para-differential calculus, we can get 

l|Gr|| B n/2—1,/2 < C ||o|| B „/2-i,„/2||divu || B „ /2 

<G(1 + t)-"/ 4 M(t)/(t), 


(5.2) 

(5.3) 

(5.4) 


l|G , 3 || B n/ 2 -l,n /2 < C ||F|| B „/2-l,n/2||Vu|| g n/2 
<G(1+t)-"/ 4 M(t)/(t), 

IG 0 1|Bn/ 2-1 <C ||aF||^„ /2 < C ||a|| B n/ 2 ||F||^„/2 

<G(1 + t)"" / 2 A/ 2 (t). 
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For the term G 2 , we need to estimates term by term carefully as follows 

(5.5) ||v • Ve|| g-n./ 2 -i + ||t; • Vt;|| An/a _ 1 < C ||«||^ /a _ 1 1|Vw||^ B/a 

(5.6) <G(l + T)- n / 4 M(T)/(T). 

Noting that (7(0) = K{ 0) = 0, we get by using Lemma 3 and Remark 6 in [5] that 

(5.7) ||G(aMu||^ /2 - 1 < G||V 2 u|| iW2 _ 1 ||G(a)|| iW2 < G(1 + r)-"%(r)/(r), 

(5.8) ||JsT(a)Va||g„ /a _ 1 < <7||Ff(a)|| i j n/2 ||Va|| i j n/2 _ 1 < G(1 + r)~ n / 2 M 2 (r), 

(5.9) II^VF'II gn/ 2-1 < C'|| J P||^„ /2 ||V J P||^„ /2 _ 1 < <7(1 + t)-"/ 2 M 2 (t). 

From the above estimates m-m, we obtain 

(5.10) || G 2 || sn /,-1 < <7(1 + t)-"/ 2 M 2 (t) + <7 (1 + t)-"/ 4 M(t)/(t). 
Substitute (l5.2D - fl5.4B and (15.101) into (15.11) . we will have 

Y E ^) < e ~° ot Y E i^ + C [ e“ Co(t_T) ((1 + t)-”/ 2 M 2 (t) 

+ (1 + r)~ n ^M (t)/(t)^ dr 

< e“ C0< Y E q(°) + M (*) / e _Co(t_r) (l + r)- n / 4 /(r)dr 

“'O 

+ M 2 (t) [ e- Co(t " T) (l + r)-"/ 2 dr 
Jo 

< e ~ Cot Y e q(°) + ^(1 + t)~ n / A 5M{t) + G(1 + t)~ n / 2 M 2 (t). 

q>R 

So we obtain 

(1 + t) n ' A Y W < C (||(a 0 , F 0 )|| B n /2 + IMIun/2-i) + CSM(t) + CM 2 (t). 

q>R 

Combining the above inequality, Remark 13.41 and using properties of homogeneous 
Besov space, we obtain 

(5.11) M(t)< C (\\(ao, Fo)\\ B n /2 + \\vo\\ B n/ 2 ~i) + C5M(t) + CM~(t). 

By taking <5 > 0 suitably small, we finally have 

(5.12) M(t) < <7(||(ao,Fo)|| iW2 + ||uo|| B „/ a -i) 

for all 0 < t < T. It follows from local well-posedness Theorem 12.81 and the above 
estimate (15.121) that 

M(t) < C < 00 

for all t > 0. Hence, we obtain the desired decay estimates in Theorem II. 21 
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